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ABSTRACT 

We study a simple model of dark matter that is gravitationally clustered around the sun in 
s^ . the form of a spherical halo of a degenerate gas of heavy neutrinos. It is shown that for neutrino 

Q , masses mv ^ 16 keV/c^, the resulting matter distribution is consistent with the constraints 

[^ ' on the mass excesses within the orbits of the outer planets, as obtained from astrometrical 

and the Pioneer 10/11 and Voyager 1/2 (Anderson et al. 1995) ranging data. However, the 
anomalous acceleration recently detected in the Pioneer 10/11 data that is approximately 
constant between 40 AU and 60 AU (Anderson et al. 1998; Turyshev et al. 1999) is incompatible 
with both our model and earlier Pioneer 10/11 ranging data for the outer planets. We then 
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\^ , calculate the planetary and asteroidal perihelion shifts generated by such a neutrino halo. For 



in 

o 



nil, ^ 16 keV/c , our results are consistent with the observational data on Mercury, Venus, 
Earth and Icarus. Finally, we propose to detect this neutrino halo directly with a dedicated 

0> . experiment on Earth, observing the X-rays emitted in the radiative decay of the heavy neutrino 

^^ ' into a light neutrino and a photon. 

^ ■ 

Q ' Subject headings: Gravitation-Dark matter-Solar system: general 
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1. Introduction 

j^ ' It is well known that the properties of galactic systems pose a great challenge to gravity theories. For 

3 , virtually all spiral galaxies the galactic rotation curves tend towards some constant value for large distances 

from the center of the galaxy. This is clearly either inconsistent with the mass distribution inferred from 
the distribution of visible stars, or, even more challenging, in contradiction with the laws of Newtonian 
dynamics. 

Most widely accepted is the conservative explanation of this discrepancy in terms of dark matter (DM) 
(e.g. Ashman 1992) which presumes that the visible stars are embedded in a massive, nearly spherical halo 
of nonluminous matter. The mass of the halo varies from one galaxy to another, but in general it constitutes 
90% of the total mass (Tremaine 1992). While the DM hypothesis could explain the flat rotation curves 
of the galaxies in a consistent manner, it has its own troubles, in particular: (i) there is no compelling 
model for the formation of the DM halo, and (ii) despite much effort, so far no known form of matter lends 
itself to a satisfactory understanding of the DM halo. For instance, the microlensing experiments MACHO 
(Alcock et al. 1996, 1997) and EROS (Ansari et al. 1996), are as yet far from explaining the galactic DM 
halo in terms of massive astrophysical compact halo objects (MACHOs). In fact, it is by no means clear 
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that the DM halo should be made of astrophysical low-mass black holes, brown dwarfs, old white dwarfs or 
more exotic nonbaryonic MACHOs. On the contrary, the DM halo may very well consist of loose clouds of 
weakly interacting (i.e. nonbaryonic) massive particles (WIMPs) such as SUSY neutralinos, neutrinos or 
axions (Jungman et al. 1996). Microlcnsing experiments have ruled out a large class of possible compact 
baryonic DM components. As gaseous baryonic components are also largely excluded, it has been argued 
that at least some nonbaryonic DM is required to explain the DM halo in our Galaxy (Freese et al. 1999). 
Of course, if DM is in the form of nonbaryonic particle clouds, it will be everywhere in our Galaxy and it 
may thus also be observed in the solar system. 

An alternative explanation for the galactic rotation curves is based on the possibility that Newton's 

gravitational law ceases to be valid at small accelerations (Milgrom 1983; Bekenstein 1992) or large scales 

(Sanders 1990). This challenging hypothesis has been discussed in a number of papers (e.g. Eckhardt 1993; 

Hammond 1994). For instance, it has been shown (Sanders 1986; Nieto and Goldman 1992) that a modified 

gravitational potential of the form 
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with a = —0.9 and rg ~ 30 kpc, could indeed explain the flat rotational curves for most of the galaxies. 
In contrast to the DM hypothesis, there would be no detectable corrections to Newton's gravity, as the 
exponential term is unity with high accuracy for the solar system. We would thus recover the standard 
form for the gravitational potential (p ~ ^GMjr within the solar system. The problem with this hypothesis 
is that it is not clear how to incorporate the modified gravitational potential in a fully general relativistic 
theory. 

The main motivation for studying DM in the solar system is that galaxies could have halos consisting of 
weakly interacting massive particles (WIMPs) ( Trimble 1987). Indeed, if DM exists in the form of WIMPs, 
one would expect this exotic form of DM to penetrate both the galactic disk and the solar system. The 
determination of the total matter density in the vicinity of the sun is a classical problem of astrophysics (see 
Tremaine (1990) and Oort (1965) for a comprehensive review of early research). It allows for a significant 
amount of DM that could have condensed into a DM halos around the stars. However, the conditions 
on particle masses and interaction cross sections under which this condensation would have taken place 
during star formation, have not been worked out in detail. In fact, it is difficult to envisage a gravitational 
interaction mechanism that would lead to capture of significant amounts of weakly interacting DM particles. 
Assuming these particles to be dissipationless, they would be part of the galactic halo and move with the 
velocities of the order of 10^"^c. In this case very few DM particles could be trapped in orbits around 
the sun and the DM density would be nearly constant near the sim. However, Bahcall (1984a, 1984b) 
has argued for the possible existence of a galactic disk component of nonluminous matter, not necessarily 
weakly interacting, with a density of the order of luminous disk matter, i.e. 0.1 Mqpc^"^. Such DM 
could have been trapped around the sun since its formation. In fact, it is possible that "nonequilibrium" 
gravitational interaction between a cloud of weakly interacting DM particles and a molecular cloud, could 
play an important role in the formation of stars with DM halos around them. The obvious place to test 
such a hypothesis is the solar system. In fact, Mikkclsen and Newman (1977) established a bound of 

Md^ 7 X IO^^Mq for DM within the Earth's orbit, while Khloper et al. (1991) made estimates of 10"'^ 
to IO^^Mq of shadow matter being captured by a baryonic matter star. If a DM halo exists around 
the sun, it would induce an additional perihelion precession of the planetary orbits (Braginsky et al. 
1992). More recently, Anderson et al. (1995) have investigated the bounds on DM within the orbits of 
different outer planets using the data from Pioneer 10/11 and Voyager 1/2. They found, within Uranus' 
orbit, an upper bound for the DM mass of Md ^ 0.5 x 10 Mq, while a boimd of Md ^ 3 x 10 Mq 



was obtained within Neptune's orbit, and finally, within Jupiter's orbit, the DM mass was determined to 
Md = (0.12 ± 0.027) X IO^^Mq. 

Thus, in this paper, we assume that massive neutrinos play the role of the weakly interacting DM 
particles, that have been trapped around the sun during its formation, and we investigate the consequences 
of the existence of such a degenerate neutrino halo. In fact, the idea that DM is made of massive neutrinos 
that can cluster around a baryonic star has been developed in a series of papers (VioUier et al. 1992, 
VioUier 1993; VioUier 1994). This model resembles the Thomas-Fermi model of an atom or ion, where the 
degenerate electron cloud around the nucleus is described by a mean electrostatic field. The difference 
between the two models is the sign of the interaction. While the model of the neutrino halo around a 
star admits also solutions without star at the center, i.e. pure neutrino balls (VioUier 1994), there is no 
corresponding solution for a nucleus with charge Z = 0. Here it is important to note that pure neutrino 
balls could also describe the supermassive compact dark objects at the galactic centers, e.g. at the center of 
our Galaxy or M87 (Tsiklauri and VioUier 1998; Bilic, Munyaneza and VioUier 1999; Munyaneza, Tsiklauri 
and VioUier 1998, 1999; Munyaneza and VioUier 1999). 

The purpose of this paper is to study the mass distribution of degenerate neutrino DM around the sun 
and to compare it with the most recent observational data from Pioneer 10/11 and Voyager 1/2 (Anderson 
et al. 1995). We will also investigate the effect of a possible neutrino halo on the precession of the perihelion 
of planetary or asteroidal orbits. 

This paper is organized as follows: In section 2, we establish the basic equations of our model and 
study the mass distribution of the neutrino halo around the sun. In section 3, we investigate the planetary 
and asteroidal perihelion shifts due to a neutrino halo. We discuss other observational implications of such 
a neutrino halo in section 4 and summarize our results with a discussion in section 5. 



2. A degenerate heavy neutrino halo around the sun 

Let us characterize the spherically symmetric cloud of degenerate heavy neutrino matter around a 
baryonic star by its gravitational potential $(r), pressure Pu{r), and mass density pu{r). In nonrelativistic 
approximation, these three quantities are linked through Poisson's equation 

A$ = 47rGp,, (1) 

the condition for hydrostatic equilibrium between the gravitational and degeneracy pressures of heavy 

neutrino matter 

dP. rf$(r) 

~T~ = ~P''^ — ' (2) 

ar ar 

and the polytropic equation of state of a degenerate nonrelativistic Fermi gas 

P. = Kpl^' . (3) 

Here the constant K is given by 
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g^ being the spin degeneracy factor of the neutrinos and antineutrinos, i.e. g^ = 2 for Majorana neutrinos 
and 5^ = 4 for Dirac neutrinos and antineutrinos. The pressure and density vanish at the radius Rq of the 
neutrino halo, i.e. 

p,(i?o) = Pu{Rq) = . (5) 



Outside the neutrino halo, r > Rq, we thus recover the standard solution of Poisson's equation (|^) 

G{M, + Mb) 
$(r) = , (6) 

where M^ is the mass of the neutrino halo, i.e. 

rRo 
M^= Airp^irydr, (7) 
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and Mb is the mass of the pointlike baryonic star at the center of the neutrino halo. 

Introducing the dimensionless potential variable v and radius x, by 



v{x)^-^^mRo)-Hr)], (8) 
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where a^ is an appropriate length scale 
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we arrive at the nonlinear Lane-Emden differential equation (VioUier, Leimgruber and Trautmann 1992, 
VioUier 1994) with index n == 3/2 

— = — fin 

subject to the boundary conditions 

"(0) = ^ and z;(xo) = 0, (12) 

Mq 

where xq — Ro/ai, is the radius of the halo in units of a^. The first condition in (03) arises from the fact 
that, near the origin, the potential energy is dominated by Mb- Thus Mb — corresponds to a pure 
neutrino ball without a pointlike source at the center. The properties of such neutrino balls have been 
discussed in a number of papers (Tsiklauri and VioUier 1998; Munyaneza,Tsiklauri and VioUier 1998,1999; 
Munyaneza and VioUier 1999). In particular, it has been shown that a neutrino ball of M « 2.6 x IO^Mq 
is a viable alternative to the supermassive black hole that many believe to exist at the Galactic center, as 
(i) it is consistent with the upper limit of the size of the supermassive compact dark object as determined 
by the motion of stars near Sgr A* (Munyaneza, Tsiklauri and VioUier 1998, 1999), and (ii) the bulk part 
of its infrared to radiowave spectrum, up to wavelengths of A = 0.3 cm, can be described by standard thin 
accretion disk theory (Bilic and VioUier 1998; Bilic, Tsiklauri and VioUier 1998; Tsiklauri and VioUier 1999; 
Munyaneza and VioUier 1999). 



The neutrino matter density may be expressed in terms of the potential variable v as 
'^^"-^ = 9.3n« [x) 



Using Eqs. (Q) and (p^), the neutrino matter mass enclosed within a radius r from the sun can be written as 

M^ir) = -Mq (xv'{x) - v{x) + v{Oy\ . (14) 

with 

M,{Ro) = -Mq ( a;ot/(xo) + i'(O) ) - M,. (15) 



Of course, the mass density p^, diverges at the origin as x~^^^ due to the pointhke nature of the 
baryonic source. However, the integral in Eq. (m) converges, yielding a finite total mass of the neutrino halo 
M^. In Fig. 1 we plot the mass A4^{r) that is enclosed within a radius r from the sun, for various values of 
the neutrino mass m,^. The slope of v at the center has been fixed to v'{0) = 1, yielding a total mass of the 
neutrino halo of M^ = 0.7 Mq. Of course, by varying f'(0), the halo could have any mass. Also shown in 

Fig. 1 is the bound on DM mass within the Earth's orbit, i.e. Md ~ 7 x 10~^Mq (Mikkelsen and Newman 
1977) as well as the constraints on the mass excesses within the orbits of various outer planets, as obtained 
from astrometrical and the Voyager 1 and 2 and Pioneer 10 and 11 ranging data (Anderson et al. 1995). 
According to these data the DM mass contained within Jupiter's orbit is Md — (0.12 ± 0.027) x 10~^Mq, 
within Uranus' orbit Md ^ 0.5 x IO'^Mq, and within Neptune's orbit Md ^ S x 10~^Mq. A DM mass 
of a few lO^'^M© over the range from 40 AU to 60 AU as shown in Fig. 1 would be consistent with the 
anomalous acceleration of ap w 7.5 x lO^^cm/s^ observed in the Pioneer 10/11, Galileo and Ulisses data 
(Anderson et al. 1998; Turyshev et al. 1999). 

Based on our model, we can determine the neutrino-mass range compatible with the limits on DM for 
the various planetary orbits. The constraints on the DM mass within Earth's orbit yields the neutrino mass 
limits 

m^c^ < 21.8 keV for g^ = 2 

m^c^ < 18.3 keV for g^ = 4. (16) 

Taking the DM data within Jupiter's orbit at face value and interpreting DM as degenerate neutrino 
matter, the neutrino mass limits are 

12.8 keV < m,,c^ < 14.2 keV for g^ = 2, 

10.8 keV < TO^c^ < 11.9 keV for g^ = 4. (17) 

However, one should perhaps interpret this range of neutrino masses as a lower limit, as Jupiter tends to 
eject any matter within its orbit (Anderson et al. 1995). 

For DM within Uranus' orbit, the upper limits for rrii, are 

m^c'^ < 12 keV for g^, = 2, 

m^c^ < 10 keV for g^ = 4, (18) 

and finally for Neptune's orbit, the bound on the neutrino mass is 

m^c^ < 16 keV for g^, — 2, 

m^c^ < 13.5 keV for g^ = 4. (19) 



We now explore what the total mass of the neutrino halo should be, in order to be consistent with the 
observational data. For this purpose, we have calculated the mass A^^(r) that is enclosed within a radius r, 
for various total masses M^ of the halo. The results are shown in Fig. 2, where the neutrino mass is fixed to 
nrii, — 14 keV/c^. As suggested by this figure, the neutrino halo cannot be heavier than ^ 10^ A/©, otherwise 
it would be inconsistent with the constraints on the observed mass excesses within the orbits of the outer 
planets. In Fig. 3, the acceleration a^, due to the presence of a neutrino halo is plotted as a function of the 
radius. Also shown is the range of allowed values for the acceleration at the positions of various planets 
(Anderson et al. 1995). We thus conclude that our model is consistent with the observations of the outer 
planets (Anderson et al. 1995). However, it cannot explain the apparent anomalous weak long-range 
acceleration seen in the Pioneer 10/11, Galileo, and Ulysses data (Anderson et al. 1998; Turyshev et al. 
1999) indicated with a horizontal bar in Fig. 3. These groups claim that there is an anomalous acceleration 
towards the sun of ap sa 7.5 x lO^^cm/s^. No variation of ap with distance from the sun was found over 
a range of 40 to 60 AU. Of course, our neutrino halo would also contribute to the anomalous acceleration 
ap. Setting ap = a^, our model predicts a decrease of the anomalous acceleration due to the fact that the 
neutrino halo mass enclosed within a distance r from the sun scales as r^^^. However, if the anomalous 
acceleration is indeed constant, as the recent observations (Anderson et al. 1998; Turyshev et al. 1999) seem 
to indicate, the enclosed mass should increase like r^ which is clearly inconsistent with the astrometrical 
and previous ranging data. Of course, there have been some attempts to explain the Pioneer anomalous 
acceleration (Katz 1998, Murphy 1998), but it is perhaps fair to say that the jury on the explanation of 
this real or spurious effect is still out. In Fig. 4, we present the mass M^ as a function of the radius Rq of 
the neutrino halo. It is interesting to note that there is always maximal radius of a halo around a baryonic 
star. For Mb = Mq, m^ = 14 keV/c^ and g^ = 2, this maximal radius turns out to be Rmax ~ 5 lyr, which 
is the scale of the interstellar distances. The corresponding neutrino halo mass would be M,j^rnax ~ 3Mq. 



3. Planetary and asteroidal perihelion shifts 

A possible neutrino halo will, of course, affect the perihelion shifts of the planets and asteroids of the 
solar system. Using Eq. (g), the gravitational potential $(r) of a sun that is immersed in a neutrino halo 
can be written as 

*(r) = $B(r)+<5$,(r), (20) 

where ^Bir) is the Newtonian potential due to the sun {Mb = Mq), i.e. 

$B(r) = -^^, (21) 

r 

and (5$^(r) is the contribution of the neutrino halo to the potential 

(5$,(r) = ^ L^ + i,'(a:o) . (22) 



Here v is the potential variable that satisfies the Lane-Emden equation with a pointlike solar mass source 
at the center, i.e. a solution of the equations (O) and (02) for Mb — Mq. 

In the vicinity of the sun, the term 5^^{r) is much smaller than <f>B(r). Near the center, the potential 
variable v{x) has the asymptotic behaviour 

v{x)Ki^x^/^ + v'{0)x + l , (23) 



where v'{0) parametrizes solutions with different halo masses. Inserting the last expression into Eq. (p2[), 
we find for the additional potential energy due to the neutrinos 

5U = m6Mr) = ll^^r^^ + C, (24) 

where C — v'{xo) — v'{0) is a constant which is irrelevant for our further considerations, and m is the mass 
of the planet or asteroid. Thus our problem reduces to investigating the perihelion shifts of the planets and 
asteroids due to the small perturbation given by Eq. (|2^) . 

When a small correction 5U{r) is added to the potential energy U = —a/r, the paths of finite motion 
are no longer closed, and at each revolution, the perihelion is displaced by a small angle (Landau and 
Lifshitz 1960) 

Here L is the angular momentum of the planet or asteroid whose unperturbed orbit is given by the standard 

equation for an ellipse, i.e. 

r = ^ , (26) 

1 + e cos if 

with 

p = a{\ — e ), e — \ A —, a = — — , a = itlGMq, L = pma, (27) 

ma^ 2£/ 

where E is the energy of the planet or asteroid (£' < 0), e is the eccentricity, a the semi-major axis and p 

the latus rectum of its orbit. 

Before calculating the perihelion shifts due to the neutrino halo, let us consider a general perturbation 
potential of the power law form 

5L/ = 7r", (28) 

where 7 and n are constant real numbers. For instance, it is well known that the general relativistic 
corrections for a planet moving around the sun can be incorporated in the Newtonian framework via an 
effective potential of the form 

^^ff = -^ + ^^^-^^ (29) 

where the constant /? is given by 

^ = ^^. (30) 

mc 

The last term in the Eq. (p9|), is the general relativistic correction to the potential SU with j — P and 
n = —3. Due to this term, the perihelion is shifted by an angle dipcRT which can be easily calculated from 
Eq. (E3) yielding 



oL 



2m f" 2-/3 



^ Jo ^ 



9 d 
= -2m^af3— 
oL 



1 r 

— J / (1 + ecos(^)d(y9 



= 6W/3aL-4 = ^^^^^ (31) 

c^a(l — e^) 

for the general relativistic perihelion shift per revolution. Let us assume that, in addition to the general 
relativistic perturbation potential ~-j3/r^ , there is a further perturbative potential given by Eq. (EN) due to 



some type of DM. It is clear that this small shift will be additive in first order. The contribution to the 
shift per revolution due to the perturbation potential SU is thus 






/ jr r dip 



= ^"^^01 
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n+2 



{1 + e COS If) 



n+2 



dip 



After some algebraic manipulations, we arrive at the expression for the shifts per revolution 



a ' 



n+l 



{2n + 3)/„+2(e) + 2^(n + 2) ( /„+2(e) - /„+3(e) 



where the function /„(e) has been defined as 

1 



In{e) = 



(l + ecos(/3)' 



-dip = 



(l-e2) 



2\"/2 



Pn- 



Vi^^ 



(32) 



(33) 



(34) 



Pn-i{x) being the Legendre polynomial of degree n — 1. In the limit of small eccentricities, i.e. e — > 0, 
Eq. (H) yields 



a 



(35) 



The perihelion shifts become negative, i.e. Sip < for ti > and n < —1 assuming of course 7 > 0. In the 
following, we will deal with two special cases: degenerate neutrino DM and DM with a constant density 
which could describe homogeneous WIMP DM in the vicinity of the sun. 

In the case of a neutrino halo, the perturbative potential SU is given by Eq. (|2j). We can thus apply 
Eq. (p3), with n = 1/2 and 7 = 4:mGMQ/{3aJ ). The constant C in Eq. (e3) does not contribute to the 
perihelion shift. After integration, we arrive at the perihelion shift per revolution 



Sip = 



8 f aV^^ (e + l)(l-e)i/2 



3 \a^ 



E{Tr/2,k) + {e- l)F{Tr/2,k) 



(36) 



where k is given by 



k^ 



2e 



1 + e 
and F(Tr/2, k) and _E(7r/2, k) are the complete elliptic integrals of the first and second kind, respectively, 

defined as 

dip 



(37) 
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V 1 — fc2 sin ip 
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E{Tr/2,k)= dipJl - fc2 sin^ ip . 
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We can now write the expression for the perihelion shift per century as 



(38) 
(39) 



Sip = -5.50" X 10^ 



E{Tr/2,k) + {e-l)F{TT/2,k) 



100 

7^' 



where Kepler's third law has been used to eliminate the period T of the planet, i.e. 

a3/2 



T = T, 



E- 



3/2 ■ 
^E 



(40) 



(41) 



Here, Te and ue are the period and the senii-niajor axis of the Earth's orbit, respectively. From Eq. (EG), 
we conclude that the shift does not depend on the semi-major axis a, i.e. they are a function of the 
eccentricity e only. Moreover, the shift is negative in contrast to the general relativistic one (see Eq. (^^). 
In Fig. 5, we present the shift due a neutrino halo as a function of the eccentricity for various neutrino 
masses. The data points show the difference between the observed perihelion shift (obs) and the correction 
due to the general relativity theory (GRT), i.e. 5(pobs — ^^grt- Thus, if there was no neutrino halo, 5(p 
should vanish identically. The general relativistic corrections, calculated using the exact general relativistic 
equations rather than the useful approximation eq. ( EQ), as well as the observed values are taken from 
Weinberg (1972). The data points with error bars are shown for Venus, Earth, Mercury and Icarus. 
As expected, the error bars of the data points decrease for increasing eccentricity. We may thus obtain 
a good upper bound for the neutrino mass using the perihelion shift data for Icarus. In Fig. 6, the 
expected neutrino halo shift is plotted for Icarus as a function of the neutrino mass. The horizontal lines 
denote the region of dtp — 5(pobs — SifcRT allowed by observations. The observational shifts for Icarus is 
Sipobs = 9.8" ± 0.8" per century, while the general relativistic correction is 5(pGRT = 10.3" per century 
(Weinberg 1972). Requiring that the shift due to a possible neutrino halo cannot be smaller than the lower 

limit of Sip = Spobs — SpGRT = —0.5" ± 0.8", i.e. Sp ^ —1.3" per century, we obtain an upper limit for the 
neutrino mass of 

TO^c^ < 16.4 keV for g„ = 2, 

m^c^ < 13.8 keV for g^ = 4 , (42) 

as seen from Fig. 6. Here, we note that the bounds (E2h agree very well with those obtained from the mass 
excesses in Figs. 1 and 2. 

Moreover, considering the violent supermassive dark object at the center of M87, with mass 
A'l — (3.2 ± 0.9) X 1O^M0 (Macchetto et al. 1997), as a relativistic neutrino ball at the Oppenheimer- 
Volkoff limit (Bilic, Munyaneza and VioUier 1999), the neutrino mass is constrained by 



m^c^ 


< 16.5 keV for .g^ 


= 2 


m^c^ 


< 13.9 keV for g^ 


= 4 



(43) 

Such a neutrino ball would be virtually indistinguishable from a supermassive black hole, as its radius is 
4.45 Schwarzschild radii, little more than the radius of the last stable orbit around a black hole of three 
Schwarzschild radii. In fact, as the mass density is very small near the surface of the neutrino ball, the 
"effective" neutrino ball radius is substantially smaller than 4.45 Schwarzschild radii. Furthermore, if we 
interpret the supermassive dark object of mass Af = (2.6 ± 0.2) x IO^A/q at the Galactic center (Chez et 
al. 1998) as a neutrino ball (VioUier, Trautmann and Tupper 1993), the observed motion of stars close to 
Sgr A* yields an upper limit for the radius of the compact dark object, and therefore a lower limit for the 
neutrino mass (Munyaneza, Tsiklauri and VioUier 1998, 1999). Finally, the infrared drop of the emission 
spectrum of Sgr A*, interpreted in terms of standard thin accretion disk theory, provides us with an upper 
limit of the neutrino mass (Bilic, Tsiklauri and VioUier 1998, Tsiklauri and VioUier 1999; Munyaneza and 
VioUier 1999), because of the cutoff of the emission of disk radiation inside the neutrino ball. Combining 
these two constraints, we obtain 

15.9 keV < TO^c2 < 25 kcV for g^ = 2, 

13.4 keV < 771^0"^ < 21 keV for g^ = 4. (44) 

Such a neutrino ball would differ substantially from a supermassive black hole, as the escape velocity from 
the center would be only about iiqo ~ 1700 km/s and the radius is about 4 x 10'' Schwarzschild radii. 
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Virtually all supermassive compact dark objects that have been observed so far at the centers of galaxies 
have masses in the range of 10^'^ to 10^'^Mq and could therefore be explained in the neutrino ball scenario. 

We now turn to the interpretation of DM in the solar system at constant density (Gr0n and Soleng 
1996) which could describe WIMP DM that is not clustered around the sun. In this case, we have n = 2 
and 7 = 2-KmGpd/S, where pd is the density of DM in the solar system. Eq. ( |33| ) yields for the shifts per 
century 

^^ = -^(l-Y/^xi^. (45) 

Here T is the Icarus period and a its semi- major axis . From Eq. ([l5|), we gather that the shifts are 
negative, as in the previous case, but they now depend on two parameters: the semi-major axis a and the 
eccentricity e. We can establish an upper limit on the density of DM using the Icarus data. Its eccentricity 
is e = 0.827, the semi- major axis a — 1.076 AU and the period T = 1.116 yr . By requiring that the DM 
shift is restricted to S(f ~ —1.3" per century, we get an upper bound on the DM density in the solar system 
of 

Pd ^ 1.5 X 10"^^ g/cm^ (46) 

Such a constant density DM model was studied by Gr0n and Soleng (1996) in the framework of general 
relativity. However, they derived an upper limit of the DM density of pd ^ 1.8 x 10~^^g/cm'^ from the 
perihelion motion of Icarus. If we use —0.8" instead of —1.3" for the lower limit of the DM shift of Icarus, 
as Gr0n and Soleng (1996) did, we obtain a density of pd ^ 9.4 x lO^^^g/cm'^ which is still larger by a factor 
of 5.2 than the upper limit obtained by Gr0n and Soleng (1996). The discrepancy between our results and 
those by Gr0n and Soleng can be traced back to the fact that those authors used for their calculations a 
perturbation expansion in terms of the eccentricity valid only for nearly circular orbits (e « 0) while Icarus 
has an eccentricity of e = 0.827. Using our value of the upper limit for pd, i.e. Eq. (|46|), one can calculate 
the total DM mass within the orbits of the various planets. Within the Earth's orbit, the DM mass is 
Md < 1.1 x 1O~^M0 which is in agreement with the bound of Md ~ 7 x 1O~*M0 obtained by Mikkelsen 
and Newman (1977). Within Uranus' orbit , we get a DM mass bound Md < 7.6 x 10~^Mo, and, finally 
for the DM mass within Neptune's orbit we have Md < 2.9 x 10~'*Mo. The last two bounds are clearly in 
conflict with the observed ephemeris, which allow only a DM mass of the order of a few times 10~^Mo, 
even within the orbit of Neptune (Anderson et al. 1995). We thus conclude that if the DM density is 
constant, the upper bound of DM within Neptune's orbit restricts the DM density in the solar system to 
Pd^ 1.5 X 10-"g/cm3. 



4. Other observational consequences of a degenerate neutrino halo 

We now turn to the question whether a neutrino halo, with properties as described in the last sections, 
can be observed in nature, using other than gravitational detection techniques. Let us consider the most 
conservative scenario, in which the Standard Model of Particle Physics (SMP) , minimally modified to 
accommodate three species of massive neutrinos, that are mixed through a leptonic Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, is basically correct at low energies. We further assume that our heavy neutrino 
is a Dirac neutrino, more specifically the mass eigenstate Vt , in the mass range between 10 keV/c^ and 
25 keV/c^, while the Vg and Vf^ are assumed to be massless. The Vr couples preferentially to the v'^ and to 
a lesser extent to the z^^ and z/' eigenstates of the charged weak interaction. In the framework of the SMP, 
the dominant decay mode of the mass eigenstate Vt- in the assumed mass range is the conventional radiative 
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decay, having a lifetime (Boehm and Vogel 1987) of 

TDi^r - i^rl) = 1.30 X 1015 p7.2keV y , ^^_ 

Here the CKM- matrix element Urvi{i — £, fJ-) with the largest modulus, presumably i — jJ., determines the 
dominant decay mode. The present limit from the CHORUS collaboration for \Uru ,\'^ is If^ri^ J^ < 3.3 x 10^^ 

for (5m^ ^ 100 eV^ (Sato 1999). One can thus safely conclude that the i^r is quasistable over the lifetime of 

the Universe. However, even though the r-neutrino is remarkably stable against radiative decay, r-neutrino 

matter is quite radioactive. It is perhaps so abundant that the X-ray flux generated by the r-neutrino decay 

rates 

1 , , 

Ui, = n^, (48) 

td 

where the neutrino (and antineutrino) number density n^ is obtained from Eq. (113) 

could be observable. In the vicinity of a pointlike baryonic star of mass Mb, i.e. v w AIb/Mq, may be 
rewritten as 

^^ mlg, ( 2GMb \"'' ..„. ,.17 ( MbRqV'^ ( rn^_\' .3 .,„. 

n,y ss T = 0.978 x 10 Qui -, cm . (50) 

6^2^H r J ^'^\MQrJ \l7.2keVj ^ ' 

Thus the number of photons with energy to^c^/2 emitted per unit time and volume is 

Although an energy of m^c^ /2 w 8 keV is equivalent to a temperature of roughly 10^ K, the X-ray flux 
from neutrino decays near the solar surface is too small to contribute significantly to nuclear synthesis in 
the sun or to maintaining the solar corona at a temperature of a few million degrees. At a distance of one 
astronomical unit (r = ue), the number density and decay rates per unit volume are 

n^ « 3.10 X 10^ V (y^^) '^^-'^ (52) 

and 

/ m c \ 
n, « -2.38 X IQ-^g, [yfi^^) \Uru.K,fcu,-^y^-\ (53) 

Thus, if our model of the solar neutrino halo is correct, we predict that in a shielded vacuum of 1000 m'^, 
one will observe g^ photons per hour for |[/i-j,^[/*^ p — 3.3 x 10"^ which is the present experimental limit of 
the CHORUS collaboration, and rrii, — 16 keV/c^, g^ being the spin degeneracy factor of neutrinos (Viollier, 
Leimgruber and Trautmann 1992). However , this number could be substantially enhanced through the 
gravitational field of the Earth. The photons originating from the the radiative decay i/,- — j i^^ + 7 or ( 
i^T- — > t^e + 7) will have a sharp energy of m,j<? /2 with l^E / E k, 10^^. This measurement could prove 
the existence of the massive neutrino halo and fix the v^ mass and the v^ — v^ ( or Vt — v^ mixing angle 
accurately. At the same time, this could be the first direct evidence for the neutrino background and the 
nature of DM. 
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If both neutrinos and antineutrinos are present in a neutrino halo around the sun, these wiU annihilate 
into light neutrinos v^, Ve and i'^, F^ through ordinary weak interactions processes via the Z", which are 
independent of the mixing angle of course. In fact, the rate of change of the neutrino (and antineutrino) 
number density is 

hy^ - < OAVy > -y, (54) 

where the spin averaged annihilation cross section is for Dirac neutrinos given by 

< (JAV„ >= — 7T^- (55) 

The largest annihilation rate is obtained in the interior of the star which of course depends on the internal 
structure of the star. However, the effective annihilation time ta at the surface of the star with mass Mb 
and radius Rb is a good indicator for how fast this process actually is. It can be calculated quite reliably 
yielding 



TA 



UiT^h' ( Rb 



< (JAVu > n^, Gtmlc \2GMb 



3/2 



. o.«...o.3(iL^)=(|^)-,. 

Thus for a neutrino mass of TO,y — 16 keV/c^ and a solar mass Mb — Mq with a radius Rb — Rq, we 
obtain ta = 0.623 x 10^'^yr, much larger than the age of the universe. Although this annihilation process 
is more efficient than the radiative decay, it will be essentially unobservable due to the low energy of the 
neutrinos. 



5. Summary 

In this paper, we have investigated the properties and implications of a possible halo of degenerate 
neutrino matter around the sun. For small halo masses or sufficiently close to the center, this neutrino halo 
is dominated by the gravitational potential of the sun. We have established that the enclosed mass of a 
degenerate neutrino halo around the sun is, for neutrino masses of TOj^ ^ 15 keV/c^, of the order of a few 
times IO^^Mq within Uranus' and Neptune's orbits, consistent with available observational data. If such 
a neutrino halo exists, it would decrease the perihelion shifts, i.e. the neutrino halo shifts are negative in 
contrast to the general relativistic ones. The perihelion shifts due to the neutrino halo depend only on one 
parameter, the eccentricity e, while those due to general relativistic effects depend on the eccentricity e and 
the semi-major axis a. The maximal radius of such a degenerate neutrino halo around the sun is a few light 
years, with a total halo mass of ^ 3 Mq. 



In order to explain the mass excesses within the orbits of various outer planets using the Voyager 1/2 
and Pioneer 10/11 data, the neutrino mass should be in a narrow range of m^ = (15 ± 1) keV/c^ for g^ = 2. 
The predicted values of the acceleration in the solar system have been compared to those obtained from 
the recent observations. We have seen that a degenerate heavy neutrino halo around the sun fits the recent 
ephemeris very well, but it cannot explain the anomalous acceleration in the Pioneer data (Anderson et al. 
1998; Turyshev 1999). We have shown that, in order to be consistent with the observations, the total mass 
of the neutrino halo should be less than ^ 1O^M0. A neutrino mass in the range from 14 to 16 keV/c^ 
fits all the observational data on DM in the solar system, the compact dark central object Sgr A* in our 
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Galaxy with M = (2.6 ± 0.2) x IO^Mq, and the most massive compact dark object at the center of the 
galaxy M87 with M = (3.2 ± 0.9) x IO^Mq. 



We have established an upper bound of DM with a constant density around the sun of 
Pd ~ 1.5 X 10~^^g/cm'^ using the upper limit for the DM mass within Neptune's orbit. 



We have proposed a new experiment aimed at observing the radiative decay of the neutrinos in the 
mass range around 15 keV/c^, which the sun might have accumulated in a degenerate neutrino halo during 
its formation. 

This research was supported by the Foundation for Fundamental Research. F. Munyaneza gratefully 
acknowledges funding from Deutscher Akademischer Austauschdienst (DAAD). We wish also to thank 
R. Lindebaum for useful discussions. 
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Figure captions: 

Figl: Mass A4^{r) of degenerate neutrino matter enclosed within a radius r from the sun for 
various neutrino masses. The total mass of the halo is 0.7Mq. The observational data points with error 
bars (Anderson et al. 1995) are shown for Jupiter (a = 5.2 AU), Uranus (a = 19.2 AU) and Neptune 
(a = 30.1 AU). The data point at r = 1 AU is taken from Mikkelsen and Newman (1977). The arrows 
indicate upper limits. The dark mass needed to explain the anomalous acceleration (Anderson et al. 1998; 
Turyshev et al. 1999) is indicated by a bar at 40 to 60 AU. 

Fig2: Mass A^^(r) of degenerate neutrino matter enclosed within a radius r from the sun for various 
total masses M^, of the neutrino halo. The neutrino mass and degeneracy factor are fixed to rrii, = 14 keV/c^ 
and g^ — 2, respectively. A total mass of the neutrino halo which is less than M^ ^ IOOMq would be 
consistent with the observed mass excess data from Pioneer 10 and 11 (Anderson et al. 1995). The data 
points with error bars arc shown for Jupiter, Uranus and Neptune. The observational data at r = lAU is 
taken from Mikkelsen and Newman (1977). The bar at 40 to 60 AU represents the dark mass needed to 
explain the anomalous acceleration of ap = 7.5 x lO^^cm/s . 

Fig3: Excess acceleration a^ due to a neutrino halo around the sun. The total mass of the halo is 
M^ — 0.7 Mq. The expected values of the acceleration using observational data points (Anderson et al. 
1995) are shown by points with arrows. The point with arrow at r = 1 AU is calculated using the data from 
Mikkelsen and Newman (1977). The anomalous acceleration at 40 to 60 AU is indicated by a horizontal 
bar. 

Fig4: The total mass M^ as a function of the radius Rq of the neutrino halo around the sun. 
The neutrino mass tti^ is varied as shown on the graph. M^^max turns out to be fairly constant, i.e. 
M^^max ~ 3M0, for the maximal radius Rmax of the neutrino halo. 

Fig5: Perihelion shifts 5(p caused by a neutrino halo as a function of the eccentricity e. The data 
points with error bars denote the difference between the observed perihelion shifts and general relativistic 
corrections for the perihelion shifts. The data points shown are for Venus (e — 0.007), Earth (e = 0.017), 
Mercury (e == 0.206) and Icarus (e = 0.827). 

Fig6: The Icarus perihelion shift Sip as a function of the neutrino mass m^j. The horizontal lines show 
the difference between the observed value and the correction predicted by general relativity theory . In 
order to be consistent with the observational data for the Icarus perihelion shift, the neutrino mass is 
constrained by 1111,0^ < 16.4 keV for g^ = 2 oy m^c^ < 13.8 keV for gi, — 4. 
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